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Abstract Spatial and/or temporal propagation of light 
waves in periodic optical structures offers a rather unique 
possibility to realize in a purely classical setting the op- 
tical analogues of a wide variety of quantum phenom- 
ena rooted in relativistic wave equations. In this work a 
brief overview of a few optical analogues of relativistic 
quantum phenomena, based on either spatial light trans- 
port in engineered photonic lattices or on temporal pulse 
propagation in Bragg grating structures, is presented. 
Examples include spatial and temporal photonic ana- 
logues of the Zitterbewegung of a relativistic electron, 
Klein tunneling, vacuum decay and pair-production, the 
Dirac oscillator, the relativistic Kronig-Penney model, 
and optical realizations of non-Hcrmitian extensions of 
relativistic wave equations. 



1 Introduction 

Quantum-classical analogies have been explored on many 
occasions to mimic at a macroscopic level many quan- 
tum phenomena which are currently inaccessible in mi- 
croscopic quantum systems [1 . The study of analogies, 
besides to be fruitful in gaining insights into different 
phenomena in nature, provides a noteworthy strategy 
in research capable of transferring ideas, concepts and 
techniques among apparently unrelated physical fields. 
In particular, in the past two decades engineered pho- 
tonic structures have provided a useful laboratory tool 
to investigate and visualize with classical optics the dy- 
namical aspects embodied in a wide variety of coher- 
ent quantum phenomena encountered in atomic, molec- 
ular, condensed-matter and matter- wave physics [2J. The 
study of quantum-optical analogies has been greatly stim- 
ulated by the development of the field of discrete op- 
tics, aimed to realize novel functionalized optical materi- 
als based on evanescently-coupled optical waveguides in 
which the dispersion and diffraction properties of light 
can be specifically managed [3 S ,4 ( 5 . The development 



of recent and reliable technologies in waveguide fabrica- 
tion, notably the one based on waveguide inscription in 
transparent glasses based on femtosecond laser writing 
[5], has enabled to access in the lab to a wide variety of 
fascinating optical analogues of quantum systems. Such 
analogies have also transferred into optics some impor- 
tant ideas and methods for molding the flow of light orig- 
inally developed in the context of quantum control [UH] . 
Among the wide variety of quantum-optical analogies 
investigated in the last two decades, we mention the op- 
tical analogues of electronic Bloch oscillations [3J[7] and 
Zener tunneling [8l[9] , dynamic localization [10] , coherent 
enhancement and destruction of tunneling [IT] , adiabatic 
stabilization of atoms in ultrastrong laser fields [12] , An- 
derson localization [T3], quantum Zeno effect [TTJ, Rabi 
flopping [IS], coherent population transfer [JJj], coher- 
ent vibrational dynamics [17], geometric potentials [18], 
and dynamical (Kapitza) trapping [JJ5]. Most of such 
analogies are based on the formal similarity between 
the paraxial optical wave equation in dielectric media 
and the single-particle nonrelativistic Schrodinger equa- 
tion, thus providing a test bed for classical analogue 
studies of nonrelativistic quantum phenomena. Optical 
analogues of many-body phenomena in nonrelativistic 
models, such as photonic analogues of the famous Bose- 
Hubbard model, have been proposed as well [20 . 

Recently, a great attention has been devoted toward 
the investigation of experimentally-accessible and con- 
trollable classical or quantum systems that simulate cer- 
tain fundamental phenomena rooted in the relativistic 
wave equations, such as the Dirac equation for fcrmionic 
particles. Among others, cold trapped atoms, ions and 
graphene have proven to provide accessible systems to 
simulate relativistic physics in the lab, and a vast liter- 
ature on this subject has appeared in the past few years 
(see, for instance, [21,22,23,51] and references therein). 
In particular, low-energy nonrelativistic two-dimensional 
electrons in graphene obey the Dirac- Weyl equation and 
behave like massless relativistic particles. This has lead 
to the predictions in condensed-matter physics and atom 
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optics of phenomena analogous to Zitterbewegung [2"5] 
and Klein tunneling [2 6) of relativistic massive or mass- 
less particles, with the first experimental evidences of 
Klein tunneling in graphene [27] , carbon nanotubes [35] , 
trapped ions [2H] > and of Zitterbewegung [30] . 

Photonic analogues of Dirac-type equations have been 
also theoretically proposed for light propagation in cer- 
tain triangular or honeycomb photonic crystals, which 
mimic conical singularity of energy bands of graphene 
[3"TII52"II55] . as well as in metamaterials [31], optical su- 
perlattices [55115511571155] . Bragg gratings [BTflHDHT] . and 
nonlinear quadratic media |42) . Such studies have moti- 
vated extended investigations on the properties of 'pho- 
tonic graphene' [5TI[55ir4"5] and to the proposals of pho- 
tonic analogues of relativistic phenomena like Zitterbe- 
wegung [321IMP5P2]. Klein tunneling [31I31IM1I22 , de- 
cay of the quantum vacuum and pair production [57] . 
the Dirac oscillator [40], and the relativistic versions of 
the Kronig-Penney model and surface Tamm states [H] . 
Noticeably, the introduction of gain and/or loss regions 
in the optical medium can be exploited to realize in a 
classical setting certain non-Hermitian relativistic mod- 
els proposed in the context of non-Hermitian quantum 
mechanics and quantum field theories [441145] . 

In this article a brief overview of a few important 
optical analogues of relativistic quantum phenomena is 
presented. Among the various optical realizations of rel- 
ativistic wave equations proposed in the recent litera- 
ture and briefly mentioned in previous discussion, the 
present work focuses mostly on a few rather simple pe- 
riodic optical structures, namely spatial light transport 
in evanescently-coupled waveguide arrays and temporal 
pulse propagation in Bragg grating structures. Specific 
examples of quantum-optical analogies discussed in the 
work include spatial and temporal photonic analogues of 
Zitterbewegung, Klein tunneling, instability of the quan- 
tum vacuum and pair-production, the Dirac oscillator, 
the relativistic Kronig-Penney model, and optical simu- 
lations of non-Hermitian relativistic wave equations. 

2 Photonic Zitterbewegung 

Originally predicted by Schrodinger in the study of the 
Dirac equation [JS], Zitterbewegung (ZB) refers to the 
trembling motion of a freely moving relativistic quan- 
tum particle that arises from the interference between 
the positive (electron) and negative (positron) energy 
states of the spinor wavefunction [251147] . For a free elec- 
tron, the Dirac equation predicts the ZB to have an 
extremely small amplitude (of the order of the Comp- 
ton wavelength ~ 10 -12 m) and an extremely high fre- 
quency (~ 10 21 Hz), making such an effect experimen- 
tally inaccessible. Moreover, the physical relevance of ZB 
in relativistic quantum mechanics is a rather controver- 
sial issue, because such an effect arises in the framework 
of the single-particle picture of the Dirac equation, but 



not in quantum field theory. Phenomena analogous to 
ZB, which underlie the same mathematical model of the 
Dirac equation, have so far predicted in a wide variety 
of quantum and even classical physical systems, as dis- 
cussed in the introduction. Here we briefly present two 
optical simulations of the relativistic ZB, based the for- 
mer on spatial light propagation in binary waveguide 
arrays, the latter on the frequency conversion process of 
optical pulses in quadratic nonlinear media. 

2.1 Photonic Zitterbewegung in binary waveguide arrays 

The simplest example in optics where one can find an 
analogue of the relativistic ZB is may be the discrete 
transport of light waves in a one-dimensional binary ar- 
ray [35] . Such an optical system behaves, in fact, like 
a one-dimensional superlattice, in which the two mini- 
bands of the superlattice plays the same role as the 
positive- and negative-energy branches of the Dirac equa- 
tion [48] . Under excitation by a broad light beam near 
the Bragg angle, the discretized light propagates along 
the array showing a characteristic trembling motion, which 
mimics the relativitsic ZB of the electronic wave func- 
tion. Such an oscillatory motion has been recently ob- 
served in Ref . [38 in binary arrays written in fused silica 
by the femtosecond laser writing technology [SJ, provid- 
ing the first simulation in optics of ZB. 

Let us consider light propagation in a binary waveg- 
uide array, realized by two interleaved lattices A and B, 
as shown in Fig. 1(a). In practice, the superlattice can be 
realized by a sequence of equally-spaced waveguides with 
alternating deep/shallow peak refractive index changes 
dni and dn2- In the tight-binding approximation, light 
transport is described by the coupled-mode equations 
[53151] 

.dc^ _ + C ;_i) + (-l) l ac h (1) 

dz 

where q are the modal field amplitudes in the vari- 
ous waveguides, 2a and n are the propagation constant 
mismatch and the coupling rate between two adjacent 
waveguides of the array, respectively. The superlattice 
supports two minibands, separated by a narrow gap of 
width 2<T [see Fig. 1(b)], defined by the dispersion curves 
io±(q) — ±y/a 2 + 4k 2 cos 2 (qa), where 2a is the lattice 
period and q the Bloch wave number. Hence, in the vicin- 
ity of the edges of the first Brillouin zone, e.g. near q = 
7r/ (2a), the dispersion curves of the two minibands form 
two opposite hyperbolas, and therefore mimic the typi- 
cal hyperbolic energy-momentum dispersion relation for 
positive-energy and negative-energy branches of a freely- 
moving relativistic massive particle (dotted graph). This 
suggests that light transport in the lattice for Bloch 
waves with wave number q close to 7r/(2a) simulates 
the temporal dynamics of the relativistic Dirac equa- 
tion. When launching a tilted broad beam E(x, z) at 
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the Bragg angle 9 b — A/(4n s a) (n s is the substrate re- 
fractive index) into the array only a small region around 
q = w /{2a) in g-space is excited. After setting C2 n {z) = 
(— \) n tpi{n, z) and c% n -\ — ~ i{~ l) n V2 z ) an d intro- 
ducing the continuous transverse coordinate £ n = 
x/{2a), the two-component spinor ip{£,z) = {tpi,tp2) T 
satisfies the one-dimensional Dirac equation [3"5H8] 



.dtp dtp . 

(— = -ik(j x — + aa z tp 
oz at; 



(2) 



where a T = 



are the a r and 



lo) and az= (o-l 
o~ z Pauli matrices. Equation (2) corresponds to the one- 
dimensional Dirac equation for a relativistic freely-moving 
particle of mass m [47: provided that the formal change 



mc 2 /h, 



(3) 



is made, and £ and z are interpreted as the spatial and 
the temporal variables, respectively. Therefore, in the 
optical superlattice the temporal evolution of the Dirac 
spinor wave function tp is mapped onto the spatial evolu- 
tion along z of the field amplitudes tpi and tp 2 , describing 
the occupation amplitudes of light in the two sub-lattices 
A and B of the binary array. Correspondingly, ZB is ob- 
served as a quivering spatial oscillatory motion of the 
beam center of mass {n){z) = ^2 n n\c n \ 2 / J2 n \°n\ 2 - Note 
that the measurable quantity (n.) (z) is directly related to 
the expectation value of position for the relativistic par- 

tide, (${z) = (/ da eci^ii 2 + i^i 2 )) / (j dam 2 + \h 2 )), 

by the simple relation (n) ~ 2(£) + 1/2 [35] ■ Let us 
indicate by E{x, 0) the electric field envelope that ex- 
cites the superlattice at the z = input. For beam 
incidence close to the Bragg angle 9b, one can write 
E{x, 0) = G{x) exp{2mxn s 9 'b / 'A) , where the envelope G 
varies slowly over the waveguide spacing a. Under such 
an assumption, an exact expression for can be de- 

rived and reads [55] 

(£)(z) = (£}{0)+v z + 2ttko- 2 [ dk{l/e 3 )sm{2ez)\G{k)\ 2 



(4) 

where k — 2qa — n is the shifted transverse momentum, 
e(fc) = \J a 2 + K 2 k 2 defines the energy-momentum dis- 
persion relation of the free relativistic particle, G{k) = 
(l/27r) J d^G{2a£) exp(— ik£) is the angular spectrum of 
the beam envelope (normalized such that J dk\G{k)\ 2 = 
l/47r), and v — Attk 3 J dk{k/e) 2 \G{k)\ 2 is the mean 
particle speed. The last oscillatory term on the right 
hand side of Eq.(4), superimposed to the straight tra- 
jectory defined by the first two terms, is the ZB. For 
G{k) centered at k = 0, at leading order Eq.(4) yields 
(£)(z) ~ (£,){0)+v o z+{k/2<t) sm(2<7z), i.e. the amplitude 
and frequency of ZB are given by 




R ZB = n/{2a) = h/{2mc) 
u>zb = 2cr = 2mc 2 /h, 



(5) 
(6) 



Fig. 1 Photonic analogue of relativistic Zitterbewegung in a 
binary waveguide array, (a) Schematic of the optical superlat- 
tice [refractive index change profile n{x) — n 3 \. (b) Dispersion 
curves of the two minibands of the superlattice, correspond- 
ing to the electron and positron energy branches of the Dirac 
equation, (c) Evolution of a broad Gaussian beam (snapshot 
of \c n {z)\ 2 , left panel), exciting the binary array and tilted at 
the Bragg angle, and corresponding beam trajectory (right 
panel), for parameter values k = 1 and a — 0.6. (d) Same as 
(c), but for a — 1. 



respectively. Therefore, ZB vanishes for either the far- 
relativistic (m — >■ 0) and the weak-relativistic (m — > 00) 
limits: in the first case the amplitude of ZB diverges, but 
the oscillation frequency luzb goes to zero, whereas in 
the latter case the frequency of ZB diverges but its am- 
plitude vanishes (see, for instance, [IE])- As an example, 
Figs. 1(c) and (d) show the ZB as obtained by numerical 
simulations of the coupled-mode equations for two val- 
ues of the detuning a and for a broad Gaussian beam 
that excites the array at the Bragg angle. The damp- 
ing of the oscillations observed in the figures is related 
to the spectral angular broadening of the incident beam 
[35j . The previous analysis can be extended to a pho- 
tonic superlattice in two transverse spatial dimensions 
x and y. The superlattice is composed by two inter- 
leaved triangular lattices, with a propagation constant 
detuning a between the guided modes, and a broad in- 
put beam is launched at the Bragg angle. In this case, 
in the continuous limit the evolution of the beam enve- 
lope along the longitudinal direction z is described by 
a two-dimensional Dirac-type equation, which reduces 
to the massless case of the previously-studied photonic 
graphene [31 , 32 , 33 in the a — » limit. 
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2.2 Photonic Zitterbewegung in nonlinear frequency 
conversion 

A second example of a photonic analogue of ZB is pro- 
vided by the process frequency conversion of optical pulses 
in a nonlinear quadratic medium arising from group ve- 
locity mismatch [42]. Let us consider the propagation 
of three optical pulses at carrier frequencies ui± , u>2 and 
LJ3 = u)\ -\- L02 in a nonlinear quadratic medium and in 
presence of group velocity mismatch. To study the ana- 
logue of ZB in the frequency conversion process, we as- 
sume that at the input plane z — the nonlinear crys- 
tal is excited by a strong and nearly continuous-wave 
pump field at frequency u\, and by a weak and short 
signal pulse at frequency L02 and temporal profile g(t). 
Assuming that group velocity dispersion is negligible 
and assuming perfect phase matching, in the no-pump- 
depletion approximation the sum-frequency generation 
process in described by the two coupled wave equations 



d_ 

dz 



j_d_ 

Vg 2 8t 



An 



d Id.. 

oz v g3 at 



—inA^ 
—inA2 



(7) 
(8) 



where At (I = 2,3) is the amplitude of the electric field 
envelope at frequency u>i, normalized such that \Ai\ 2 (I — 
2,3) is the photon flux at frequency u>i, v g i is the group 
velocity in the medium at frequency u)i, k = p\J Ii/Tllui, 
p is the strength of the nonlinear interaction, and I\ is 
the intensity of the pump field. The coupled equations 
(7) and (8) can be cast in a Dirac form after introduction 
of the coordinates of a moving frame 

£ = z , rj = t- z/v g (9) 

where the velocity v g is defined by the relation 



11/1 1 

v g 2 \v g2 v g3 



(10) 



In the moving frame, after introduction of the spinor 
wave field tp = (A2, A^) T Eqs.(7) and (8) can be cast in 
the form of a Dirac equation 



.dip . dip 
1— = -ia z b— + Ka x ip 
a£ on 



(11) 



where o~ x and o~ z are the Pauli matrices and where we 
have set 

2 W V g3 J 

Note that, after the formal change 
S -> c 



(13) 




propagation distance z (mm 



z (mm) 



^ — > t , Tj — > X 



Fig. 2 Zitterbewegung in the process of sum frequency gen- 
eration induced by group velocity mismatch, (a) and (b) show 
the evolution of the signal and sum-frequency pulse intensi- 
ties in a 15-mm-long PPLN crystal for a continuous-wave 
pump. The crystal is excited by a Gaussian signal pulse of 
duration t p = 0.5 ps. In (c) the evolution of the pulse center 
of mass (77)1 for the signal wave is depicted. The oscillatory 
behavior of {n)i, arising from the group velocity mismatch 
between signal and sum-frequency fields, is a signature of ZB. 



Eq.(ll) corresponds to the one-dimensional Dirac equa- 
tion for a relativistic particle of mass m in absence of 
external fields, moving alon the x axis, written in the 
Weyl representation |47) . Such a representation is dif- 
ferent, though equivalent, to the one used in the previ- 
ous subsection [compare Eq.(ll) with Eq.(2)], the spinor 
components in the two representations being related by 
a unitary transformation. Therefore, the temporal evolu- 
tion of the spinor wave function tp for the Dirac particle 
is mapped into the spatial evolution of the envelopes A2 
and A3 for signal and sum- frequency pulses, respectively, 
whereas the spatial coordinate of the Dirac particle is 
mapped into the retarded time 77 of the optical pulses. 
In the absence of group-velocity mismatch, the mass- 
less Dirac equation is obtained, in which ZB vanishes. 
As an example, Fig. 2 shows ZB in the process of sum 
frequency generation which applies to a periodically- 
poled Liyhoum Niobate (PPLN) crystal with a quasi- 
phasc matched grating that ensures phase matching at 
the wavelengths Ai = 1550 nm (pump), A2 = 810 run 
(signal), and A3 = 532 nm (sum- frequency) . The crystal 
is pumped by a near continuous-wave pump of inten- 
sity I — 1 MW/cm 2 and excited at the input plane by 
a short Gaussian signal pulse of duration t p — 0.5 ps. 
The mismatch of the group velocities of signal and sum- 
frequency fields clearly results in an oscillating motion 
(jitter) of the signal pulse along the propagation in the 
nonlinear crystal, which is the signature of ZB. 



3 Photonic analogues of Klein tunneling 

A remarkable prediction of the Dirac equation is that 
a below-barrier electron can pass a large repulsive and 
sharp potential step without the exponential damping 
expected for a non-relativistic particle. Such a trans- 
parency effect, originally predicted by Klein |26j and re- 
ferred to as Klein tunneling (KT), arises from the exis- 
tence of negative-energy solutions of the Dirac equation 
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and requires a potential step height AV oi the order 
of twice the rest energy mc 2 of the electron [49 . Rela- 
tivistic tunneling across a smooth potential step, which 
describes the more physical situation of a constant elec- 
tric field E in a finite region of space of length I, was 
subsequently studied by Sauter [50] ■ Sauter showed that 
to observe barrier transparency the potential increase 
AV ~ eEl should occur over a distance / of the order or 
smaller than the Compton wavelength Ac =%/ {mc), the 
transmission probability rapidly decaying toward zero 
for a smoother potential increase [49,50,51 . The re- 
quired field corresponds to the critical field for e + e~ 
pair production in vacuum, and its value is extremely 
strong making the observation of relativistic KT for elec- 
trons very challenging. Therefore, growing efforts have 
been devoted to find experimentally accessible systems 
to investigate analogs of relativistic KT, as discussed 
in the introduction section. In particular, great interest 
has suscitated the proposal and first experimental evi- 
dences of KT for non-relativistic electrons in graphene, 
which behave like massless Dirac fermions. In optics, 
several proposals of KT analogs have been suggested as 
well in the past few years, as discussed in the introduc- 
tion section. An important optical set up is provided by 
light transport in honeycomb photonic lattices, the co- 
called photonic graphene, whose band structure mimics 
the one of graphene |33| . However, even in simpler one- 
dimensional periodic optical structures one can realize 
optical analogues of KT [361I39] . Here we briefly discuss 
two of such simple realizations. The former one is based 
on monochromatic light propagation in waveguide super- 
lattices |36] , whereas the latter one is based on temporal 
pulse propagation in fiber Bragg gratings [35] ■ 



3.1 Photonic Klein tunneling in optical superlattices 

As in Sec. 2.1, let us consider the transport of monochro- 
matic discretized light waves in a binary array of waveg- 
uides, realized by two interleaved lattices A and B. As 
compared to the superlattice used to realize the ZB [Fig. 1(a)] 
here a weak modulation Ri of the index change, much 
smaller than either dn\ and dn2 and equal for the inter- 
leaved lattices A and B, i.e. such that R21-1 = R21, is 
superimposed to the lattice, as shown in Fig. 3(a). Light 
transport in the superlattice is then described by the 
coupled-mode equations [compare with Eqs.(l)] 



(a) 



.da 



-k(q+i + ci-i) + (-I)'o-q + (14) 



where the weak modulation Ri of the refractive index 
change is responsible, at leading order, to a slight change 
<Pi of the modal propagation constants in the various 
waveguides. For a broad beam exciting the superlattice 
tilted close to the Bragg angle, in the continuous limit 
light transport turns out to be described by the Dirac 



ABABABABAB 

tsiiifl 

a a c u c, c l+ , 




Fig. 3 Photonic analogue of relativistic Klein tunneling in 
a binary waveguide array with an interface, (a) Schematic 
of the optical superlattice. (b) Beam refraction across a po- 
tential step, corresponding to KT. (c) Absence of beam re- 
fraction, corresponding to inhibition of KT, for a smooth 
potential barrier. The right panels in (b) and (c) depict the 
space-energy band diagram, clearly explaining the inhibition 
of tunneling in (c). 



wave equation [compare with Eq.(2)] 



dz 



-wu x ^- + aa^ + » ( 15 ) 



where = $21 = $21-1- Note that Eq.(15) is formally 
analogous to the one-dimensional Dirac equation for an 
electron of mass m in presence of an electrostatic po- 
tential $(x), once the formal substitutions (3) are made 
(see, for instance, [47). Therefore, by engineering the 
refractive index depths of the various waveguides in the 
array, a step-like potential <?(£) at £ = can be realized. 
KT can be observed as a beam transmission (refraction) 
across the interface [36] . By varying the steepness of the 
potential step, the optical superlattice can be designed 
to demonstrate the disappearance of KT in a smooth 
potential step as well, i.e. the transition from KT for a 
steep potential to the Sauter's tunneling suppression for 
a smooth potential step. In the optical context, Sauter's 
inhibition of light tunneling at a smooth interface can 
be explained by considering the space-dependent band 
structures of the tight-binding superlattices, shown in 
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Fig. 3. Such bands are similar to the energy band dia- 
grams of a semiconductor, in which the two superlat- 
tice minibands (i.e. the electron and positron energy 
branches of the Dirac equation) play the role of the con- 
duction and valence bands and the effect of the external 
potential is to curve the band structure. For a sharp 
potential step [Fig. 3(b)] at the interface £ = the elec- 
tron and positron energy bands are overlapped, and the 
beam does not need to cross any forbidden region, i.e. 
KT occurs. However, owing to the sharp discontinuity 
of media properties, beam transmission is not complete, 
and some light is reflected like in a Fresnel discontinu- 
ity between two different dielectric media. Conversely, 
for a smooth potential step [see Fig. 3(c)] the beam has 
to cross a forbidden region, which behaves like a poten- 
tial barrier. The width of the potential barrier is indi- 
cated by the segment AB in Fig. 3(c). Because the width 
AB increases as the potential step gets smoother, the 
corresponding tunneling probability, i.e. beam transmis- 
sion across the interface, rapidly decreases. This explains 
the inhibition of KT and the transition to the Sauter's 
regime [50] . 



3. 2 Photonic Klein tunneling in fiber Bragg gratings 

Another rather simple optical realization of Klein tun- 
neling is provided by photonic tunneling in structured 
fiber Bragg gratings (FBGs) 39 . In a FBG, the effective 
index of the fiber is modulated along the longitudinal z' 
direction according to n(z') — no+An m(z') cos[2irz' /A+ 
2<p(z')], where no is the effective mode index in absence 
of the grating, An <C no is the peak index change of the 
grating, A is the nominal period of the grating defining 
the reference frequency ujb — ttc/ (An ) of Bragg scatter- 
ing, c is the speed of light in vacuum, and m(z'), 2<j>(z') 
describe the slow variation, as compared to the scale 
of A, of normalized amplitude and phase, respectively, 
of the index modulation. Note that the local spatial fre- 
quency of the grating is k(z') = 2n/ A+2(d(f>/dz'), so that 
the local chirp rate is C = dk/dz' = 2(d 2 (f)/ dz' 2 ). The 
periodic index modulation leads to Bragg scattering be- 
tween two counterpropagating waves at frequencies close 
to ujb- By letting E{z' , t) = <pi(z' , t) exp[—iuiBt+ikBz' + 
i<p(z')]+if2(z' , t) exp[—iujBt — ikBz' — i4>(z')]+c.c. for the 
electric field in the fiber, where ks = n/A, the envelopes 
ipi and if2 of counterpropagating waves satisfy coupled- 
mode equations. After introduction of the dimensionless 
variables z — z' jZ and t = t/T , with characteristic 
spatial and time scales Z = 2ng/ (ksAn) and T = Z/v g , 
and the new envelopes ipi^iz') = [fi(z') =p f2(z')]/V2, 
the coupled-mode equations can be cast in the Dirac 
form 139 



id T t/j = —ia x d z ip + m(z)a z ip + V(z)ip (16) 

for the spinor wave function ip — (ipi,ip2) T , where 
V(z) = (d(j)/dz), v g is the group velocity at the Bragg 
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Fig. 4 Photonic realization of relativistic Klein tunneling in 
a structured FBG. (a) Schematic of the local period detun- 
ing V = (d(f)/dz) (upper plot) and amplitude m (lower plot) 
profiles of a FBG to observe the optical analogue of KT. 
(b) and (c): Space-energy band diagrams of the FBG for a 
sharp [(b)], and a smooth [(c)] potential step V(z) of height 
AV . The shaded regions are the forbidden (stop band) en- 
ergies, the dotted horizontal line is the frequency fio of the 
incoming wave packet, and the dashed curve is the shape 
of the potential step V(z). (d) Numerically computed spec- 
tral transmittance of three FBGs for increasing length of the 
chirped region (from left to right). The transmission window 
inside the stop band of the grating for a sharp potential step 
(left panel) is the signature of KT. 



frequency, and a XiZ are the Pauli matrices. In its present 
form, Eq.(16) is formally analogous to the one-dimensional 
Dirac equation with % = c = 1 in presence of an external 
electrostatic potential V(z), m playing the role of a di- 
mensionless (and generally space-dependent) rest mass 
(see, for instance, |491l5T] b The optical analogue of the 
forbidden energy region of the Dirac equation, in this 
case, simply corresponds to the photonic stop band of 
the periodic grating. As the refractive index modulation 
of the grating, i.e. the mass term m in the Dirac equa- 
tion (16), is decreased, the stop band region shrinks and 
the limit of a massless Dirac equation (similar to the one 
describing the dynamics of electrons in graphene near a 
Dirac point) is attained. The additional external poten- 
tial V in Eq.(16), related to the phase modulation of the 
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grating according to V(z) = (d(f>/dz), changes the local 
position of the forbidden energy region. Therefore, pulse 
propagation in a FBG with a suitably designed phase 
and amplitude grating profiles can be used to mimic the 
relativistic tunneling of a wave packet in a potential step 
V(z). A typical diagram of a structured FBG suited to 
realize KT and the inhibition of KT for a smooth poten- 
tial step, together with the space-energy diagram of the 
grating [33], is shown in Figs.4(a-c). The signatures of 
KT and of its inhibition for a smooth potential step are 
simply revealed by spectrally-resolved transmission mea- 
surements of the grating. In fact, for a sharp potential 
step a transmission window is opened inside the bandgap 
of the periodic grating, whereas for a smooth potential 
step tunneling is prevented [according to the diagram 
of Fig. 4(c)] and the transmission window inside the gap 
vanishes. This is shown in Fig. 4(d), which depicts typical 
transmission spectra in three FBGs with different chirp 
rates, realized to operate near the A = 1.5 (im wave- 
length of optical communications (for more details see 

EU). 



4 Optical simulation of vacuum decay and pair 
production 

Electron-positron pair production due to the instabil- 
ity of the quantum electrodynamics (QED) vacuum in 
an external electric field is another remarkable predic- 
tion of Dirac theory and regarded as one of the most 
intriguing nonlinear phenomena in QED, whose experi- 
mental observation is still lacking (see, for instance, [521 
153)1. In intuitive terms and in the framework of one- 
particle Dirac theory, the pair production process can 
be simply viewed as the transition of an electron of the 
Dirac sea occupying a negative-energy state into a fi- 
nal positive-energy state, leaving a vacancy (positron) 
in the negative-energy state. There are basically two 
distinct transition mechanisms: the Schwinger mecha- 
nism, induced by an ultrastrong static electric field, and 
dynamic pair creation, induced by time-varying elec- 
tric fields. The Schwinger mechanism [54) is basically 
a tunneling process through a classically forbidden re- 
gion, bearing a close connection to Klein tunneling dis- 
cussed in the previous section. The other mechanism, 
namely dynamic pair creation, was originally proposed 
by Brezin and Itzykson [55 for oscillating spatially ho- 
mogeneous fields, and has recently attracted great at- 
tention because of its potential implementation using 
counter-propagating ultrastrong laser fields. Dynamic pair 
production is closely related to such intriguing phenom- 
ena like Rabi oscillations of the Dirac sea (see, for in- 
stance, [S3]). In the framework of the one- particle Dirac 
theory of vacuum decay, a simple picture of pair produc- 
tion is represented by the time evolution of an initially 
negative-energy Gaussian wave packet, representing an 
electron in the Dirac sea, under the influence of an os- 



cillating electric field [55] . When the e + e~ pair is pro- 
duced, a droplet is separated from the wave packet and 
moves opposite to the initial one |56) . The droplet is a 
positive-energy state and represents the created electron. 
In a recent work 37 , it was shown that light trans- 
port in a binary waveguide array with a sinusoidally- 
curved optical axis provides a classical simulator of the 
dynamic pair production process, in which pair produc- 
tion is visualized as a breakup of an initial Gaussian wave 
packet, composed by a superposition of Bloch modes of 
the lowest lattice miniband and representing an electron 
in the Dirac sea. Periodic axis bending of the waveg- 
uides mimics the effect of an external ac field, which 
induces transitions into the upper lattice miniband (the 
electron energy branch). The optical structure proposed 
to simulate dynamic pair production in basically the bi- 
nary superlattice of Fig. 1(a), however the optical axis of 
the waveguides is now periodically bent along the prop- 
agation direction z. If the array is excited by a broad 
beam tilted close to the Bragg angle, the continuous 
limit of the coupled-mode equations now leads to the 
following Dirac equation for the spinor ip [37] [compare 
with Eq.(2)] 



.dtp 
dz 



dtp 



iKcr x -—- + aa z tp — 2K<P(z)a x tp, 



(17) 



where <&(z) accounts for the waveguide axis bending and 
reads explicitly [37] 



*(*) = 



2irn s a(dxo/dz) 
A 



(18) 



In the previous equation, a is the distance between two 
adjacent waveguides, n s is the bulk refractive index, A 
the wavelength of light, and xo(z) is the axis bending 
profile. Note that, after the formal change 



K 

a 

z 



mc 
~ 

x 

eA x 

2hc 
t, 



(19) 



Eq.(17) corresponds to the one-dimensional Dirac equa- 
tion for an electron of mass m and charge e in presence 
of a spatially-homogeneous and time-varying vectorial 
potential A = (A x , 0, 0), which describes the interaction 
of the electron with an external oscillating electric field 
E x (t) = —(dA x /dt) in the dipole approximation (see, 
for instance, [47] ) . Because of momentum conservation, 
in the spatially-homogeneous and time-dependent field 
the problem of pair creation can be reduced to the tran- 
sition between two states consisting of a negative and 
a positive energy state coupled by the external field. 
Within the two-level model, pair production generally 
occurs as a multiphoton resonance process enforced by 
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Fig. 5 Photonic realization of dynamic e e~ pair produc- 
tion in a binary waveguide array with bent optical axis. The 
figure shows the probability of pair production after the ap- 
plication of a single-cycle ultrastrong laser pulse for increas- 
ing values of the amplitude ^o- The insets show the detailed 
beam evolution along the array corresponding to the ampli- 
tudes of A, B and C. The fractional beams deflected at oppo- 
site angles correspond to states in the negative (positron) and 
positive (electron) energy branches of the Dirac equation. 

energy conservation, with interesting effects such as Rabi 
oscillations of the quantum vacuum. Pair production in- 
duced by an ultrastrong and ultrashort laser pulse can be 
mimicked by assuming, for instance, a single-cycle of the 
ac field [<P{z) = $ sin(27rz/yl) for < z < A, <P(z) = 
for z < and z > A]. As an example, Fig. 5 shows the 
behavior of the transition probability P, after the inter- 
action with a single-cycle pulse, as obtained by numerical 
simulations of the Dirac equation for k — 2, a = 1.817, 
A = 0.6676, momentum q = n/(4a) ,and for increas- 
ing values of the the field amplitude <Pq. The detailed 
beam evolution of an initial Gaussian wave packet along 
a waveguide array with single-cycle modulated axis, cor- 
responding to the three conditions A,B and C of Fig. 5, 
are depicted in the insets of the figure. The input wave 
packet is mostly composed by Bloch modes belonging 
to the lowest miniband of the array, i.e. it belongs to 
the Dirac sea (the negative-energy branch of the Dirac 
equation). Owing to the modulation, excitation of the 
upper moninband, corresponding to the creation of a 
e + e~ pair, is clearly observed. Note that, as the wave 
packets belonging to the two minibands refract at differ- 
ent angles, they separate each other after some propaga- 
tion distance. Such a splitting is precisely the signature 
of e + e~ pair production, as discussed in |56j . 

5 Optical simulation of the Dirac oscillator 

The relativistic extension of the quantum harmonic os- 
cillator, the so-called Dirac oscillator (DO) 57,58,59 , 
provides a paradigmatic and exactly solvable model in 
relativistic quantum mechanics. Originally proposed in 



quantum chromodynamics in connection to quark con- 
finement models in mesons and baryons |60j . the DO has 
received great interest in relativistic many-body theories 
and supersymmctric relativistic quantum mechanics (see 
[55115*^110 155] and references therein). The DO model 
is obtained from the free Dirac equation by the introduc- 
tion of the external potential via a non-minimal coupling 
[5711591 63 . Since the resulting equation is linear in both 
momentum and position operators, in the nonrelativis- 
tic limit a Schrodinger equation with a quadratic poten- 
tial is then obtained. In spite of the great amount of 
theoretical studies, the DO model in relativistic quan- 
tum mechanics and particle physics remains far from 
any experimental consideration. In a recent work [40j , a 
photonic realization of the DO, based on light propaga- 
tion in structured FBGs, has been proposed. The main 
idea is that, as shown in Sec. 3. 2, coupled-mode equa- 
tions describing forward and backward light coupling in 
a Bragg grating structure has the form of a Dirac equa- 
tion. Using the same notations as in Ref. [liU] , let us con- 
sider light propagation in a FBG with an index grating 
n(z) — hq + An h{z) cos[27rz//l + 4>{z)], where no is the 
effective mode index in absence of the grating, An <C no 
is a reference value of the index change of the grating, 
A is the nominal grating period defining the Bragg fre- 
quency uib — nc/(Ano), c is the speed of light in vac- 
uum, and h(z), 4>(z) describe the amplitude and phase 
profiles, respectively, of the grating. To study Bragg scat- 
tering of counterpropagating waves at frequencies close 
to ujbi let E(z,t) = {u(z, t) exp(— icjgi + 27rzzno/Ao) + 
v{z, t) exp(— icost — 27rzzno/Ao) + c.c.} be the electric 
field in the fiber, where Ao = 2ngA is the Bragg wave- 
length, and the envelopes u and v of counterpropagat- 
ing waves satisfy coupled-mode equations. After intro- 
duction of the dimensionless variables x = zjZ and 
r = t/T , with characteristic spatial and time scales 
Z = Xo/inAn) and T = Z/v g (v g is the group ve- 
locity at the Bragg frequency), and the new envelopes 
"01, 2(2) = [ u (z) T v(z)]/\/2, the coupled-mode equations 
can be cast in the Dirac form [40] 

id T ip = a x {p x - if{x)a z } tp + a z m(x)ip (20) 

for the spinor wave function ip = (ipi,ip2) T , where a x 
and cr z are the Pauli matrices, p x = —i(d/dx), and 

m(x) = h(x) cos[cf>(x)] , f(x) = — h(x) sm[</>(x)]. (21) 

In its present form, Eq.(20) is analogous to the one- 
dimensional Dirac equation 1471 , written in atomic units 
(h = c = 1), with a space dependence mass m and 
with the momentum operator p x substituted with p x — 
if(x)a z . The space dependence of the particle mass m is 
known to describe the particle interaction with a scalar 
Lorentz potential, whereas the substitution p x — » p x — 
if(x)a z corresponds to a non-minimal coupling which 
is essential to describe the relativistic DO (see, for in- 
stance, [S3])- To realize the one-dimensional analog of 
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the DO, let us choose the amplitude h and phase ip pro- 
files of the grating such that h cos <j> = mo and h sin </> = 
—f(x) = u> s mox, i.e. [see Fig. 6(b)] 



h(x) = moy/l + {<+> s x) 2 , <p{x) = atan(o; s a;), (22) 

where mo and uj s are two arbitrary constants, corre- 
sponding to the particle rest mass and oscillation fre- 
quency of the DO in the non- relativistic limit [63]. An- 
alytical expressions of the energy spectrum for the one- 
dimensional DO and of corresponding bound states can 
be derived following a standard procedure detailed e.g. 
in Ref . |63j . Assuming for the sake of definiteness uj s > 0, 
the positive-energy spectrum (electron branch) of the 
DO is given by 

S n = yjml + 2m w s (l + n) , n = 0, 1, 2, 3, ... (23) 

whereas the negative-energy spectrum (positron branch) 
is given by 



ml + 2moui s n , n = 0, 1, 2, 3, .. 



(24) 



The corresponding eigenfunctions ip±(x) can be sim- 
ply expressed in terms of Hermite polynomials multi- 
plied by a Gaussian function. Note that the negative 
(positron) energy spectrum of the DO is not obtained 
from the positive (electron) energy spectrum by sign re- 
versal (S — > —5), the positron branch possessing an addi- 
tional bound state with energy S = —mo. In our FBG re- 
alization of the DO, bound-states with positive and neg- 
ative energies should correspond to trapped light states 
in the FBG with resonance frequencies above (5 > 0) 
and below (<5 < 0) the Bragg frequency wg, respectively. 
In a FBG of finite length, the ideal amplitude profile 
h(x), defined by Eq.(22), must be truncated, i.e. one has 
h(x) = for |x| > L/(2Z), where L is the grating length 
[see Fig. 6(b)]. The effect of grating truncation is twofold. 
First, the bound states of the DO become actually reso- 
nance modes with a finite lifetime, which should be thus 
observable as narrow transmission peaks embedded in 
the stop band of the grating. Second, the number of res- 
onance modes sustained by the grating is finite owing 
to grating truncation. As an example, Fig. 6(c) shows a 
typical transmission spectrum [power transmission ver- 
sus normalized frequency detuning 6 = (lj — ujb)T], for a 
FBG with length L/Z = 10 for parameter values mo = 1 
and w s = 0.5. In the figure, the corresponding reflection- 
band diagram (dashed area) of the grating in the (x, S) 
plane, together with the electron (S > 0) and positron 
(5 < 0) levels and a few eigenstates of the DO, are also 
depicted. The transmission peaks visible in the spectrum 
of Fig. 6(c), embedded in the stop band of the FBG, 
occur precisely at the values S n predicted by Eqs.(23) 
and (24) . Note the asymmetry of the transmission spec- 
trum around <5 = 0, with an additional resonance in the 
positron (6 < 0) branch of the spectrum with no coun- 
terpart in the electron (d > 0) branch. In an experi- 
ment, the resonant states of the DO and the asymmetry 
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Fig. 6 Photonic realization of the Dirac oscillator in a fiber 
Bragg grating, (a) Schematic of a FBG. (b) Example of 
FBG amplitude and phase profiles that realize the analog 
of the Dirac oscillator (mo = 1, u> s = 0.5, L/Z = 10). (c) 
Numerically-computed power transmission (dB units) of a 
lossless FBG with amplitude and phase profiles shown in 
(b) (left panel), and corresponding reflection band diagram 
[dashed area, right panel] with a few low-order intensity pro- 
files of trapped modes in the electron and positron branches. 



of the spectrum can be simply detected from spectrally- 
resolved transmission measurements of the grating using 
standard techniques. 



6 Photonic realization of the relativistic 
Kronig-Penney model and relativistic surface 
Tamm states 

The Kronig-Penney model for the non-relativistic Schrodinger 
equation |64j is one among the simplest models in solid- 
state physics that describes the electronic band structure 
of an idealized one-dimensional crystal. Relativistic ex- 
tensions of the Kronig-Penney model (also referred to 
as the Dirac-Kronig-Penney model) have been discussed 
by several authors (see, for instance, [65,66 67,68, 69 , 70 , 
[7T1[7^[7^ [74 l [75 l [7o 1 and references therein), and the im- 
pact of relativity on the band structure and localization, 
such as shrinkage of the bulk bands with increasing band 
number, have been highlighted on many occasions. In 
earlier studies, the Dirac-Kronig-Penney model also at- 
tracted some attention and caused a lively debate about 
the existence of so-called Dirac surface states, i.e. rela- 
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tivistic surface Tamm states which disappear in the non- 
relativistic limit pl[70l[77l[78ll79ll80ll8Tl[52] . In Ref.[4T]. 
a photonic realization of the relativistic Kronig-Penney 
model and relativistic surface Tamm states has been 
proposed, which is based on light propagation in super- 
structure FBGs with phase defects. Light propagation 
in a FBG with a periodic sequence of phase slips was 
shown to simulate the relativistic Kronig-Penney model, 
the band structure of which being mapped into the spec- 
tral transmission of the FBG. Similarly, a semi-infinite 
FBG with phase defects interfaced with a uniform FBG 
with a different modulation period was shown to sup- 
port Tamm surface states analogous to the relativistic 
Tamm states. Such surface states are responsible for nar- 
row resonance peaks in the transmission spectrum of the 
grating. 

Using the same notations as in Sec. 3. 2, light propa- 
gation in a superstructur FBG is described by the Dirac 
equation (16), where the scalar (m) and vectorial (V) 
potential terms entering in the equation are defined by 
the apodization and phase profiles of the grating, respec- 
tively. The Dirac-Kronig- Penney model for an infinitely- 
extended lattice corresponds to a constant mass m(z) = 
mo and to a potential V(z) given by the superposition 
of equally-spaced <5-like barriers, namely [65] 



V(z) = V a 



E 

n= — oo 



S(z — net) 



(25) 



where Vq > is the area of the barrier and a is the lattice 
period. Stationary solutions i/>(z,r) = ipo( z ) ex P(~ iEr) 
to the Dirac equation (16) in the periodic potential (25) 
with energy E are of Bloch-Floquet type, i.e. ipo(z + a) = 
ipo{z) exp(iqa), where q is the Bloch wave number which 
varies in the first Brilloiun zone (— ir/a < q < ir/a). The 
corresponding energy spectrum is composed by a set of 
allowed energy bands E = E(q), which are defined by 
the following implicit equation (see, for instance, [80] ) 

cos(ga) = cos(Vo) cos(ko) H sin(Vo) sin(Ka), (26) 

K 



where we have set 



E 2 



(27) 



Equation (26) defines the dispersion relation of the rel- 
ativistic Kronig-Penney model, which has been investi- 
gated by several authors (see, for instance, 65.671170]). 
The ordinary non-relativistic limit of the Kronig-Penney 
model is attained from Eqs.(26) and (27) for Vo -C 1 
and for energies E close the mo, for which the energy- 
momentum relation (27) reduces to the non-relativistic 
one [E ~ mo + k 2 /(2too)]; in this regime, the dispersion 
relation (26) reduces to cos(qa) = cos(na)-\-(moVo/ n) sm(na 
which is the ordinary dispersion relation encountered in 
the non-relativistic Kronig-Penney model. For larger en- 
ergies E but still for a low barrier area Vo < 1, non- 
relativistic effects come into play as perturbativc effects, 



which modify positions and widths of the allowed en- 
ergy bands. Non-relativistic effects deeply modify the 
band structure of the crystal for potential strengths Vq 
of the order ~ 1. In particular, if Vq is an integer mul- 
tiple of 7r, all band gaps disappear and the dispersion 
relation reduces to the one of a relativistic free particle 
[Eq.(16) with V(z) = 0], as if the ^-barriers were absent. 
In a FBG, the Dirac-Kronig-Penney model simply cor- 
responds to an infinitely-extended uniform FBG with a 
superimposed periodic sequence of lumped phase slips 
of equal amplitude A<j) = 2Vq and spaced by the dis- 
tance a. The circumstance that the effects of the 8 bar- 
riers disappear in the Dirac-Kronig-Penney model when 
Vq is an integer multiple tt is simply due to the fact 
that, under such a condition, the phase slips are integer 
multiplies of 2tt, and thus the grating has no phase de- 
fects and mimics the dynamics of a one-dimensional free 
relativistic Dirac particle. A special case of the Dirac 
Kronig Penney model corresponds to case Vo = tt /2. In 
this case, the resulting superstucture FBG comprises a 
periodic sequence of tt phase slips, which has been pro- 
posed and demonstrated to realize slowing down of opti- 
cal pulses [83] . The band structure of the Dirac-Kronig- 
Penney model is simply mapped into the alternation of 
stop/transmission bands observed in spectrally- resolved 
transmission measurements of the FBG. An example of 
the spectral transmission features of a superstructure 
FBG with periodic tt phase slips, in which the allowed 
spectral transmission bands reproduce the band struc- 
ture of the Dirac Kronig Penney model, is shown in 
Fig.7. 

If the periodic potential V(z) is truncated, surface 
Tamm states do appear. Such states for the relativis- 
tic Kronig-Penney model attracted some interest in ear- 
lier papers by several authors [6^1l70ll77ll78l[79ll80ll8Tll82] . 
and a lively debate was raised about the proper bound- 
ary conditions that should be imposed to the relativistic 
wave function at a 8 barrier. As earlier works [5511771 
1751175] suggested that the relativistic treatment yields 
a new class of surface states (the so-called Dirac sur- 
face states) which do not correspond the common Tamm 
states in the non-relativistic limit, it was subsequently 
realized that application of more physical boundary con- 
ditions does not yield any surface state which violates 
the Tamm condition in the non-relativistic limit (80) . 
The relativistic extension of the Tamm model is defined 
by the potential (see, for instance, [701150] ) 



V(z) 



Vi 



z < 
na) z > 



(28) 



Surface states are found as localized solutions to Eq.(16), 
near the surface z = 0, satisfying the appropriate bound- 
ary conditions, as discussed in [80] . In our photonic sys- 
tem, the potential V(z) defined by Eq.(28) and support- 
ing the surface Tamm states at the z = boundary is 
basically realized by two adjacent sections of uniform 
grating regions but with different grating periods, with 
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Fig. 7 Photonic realization of the Dirac-Kronig-Penney 
model in a superstructure FBG comprising a periodic se- 
quence of 7T phase slips. The figure shows the spectral power 
transmission of the grating with amplitude and phase profiles 
given in the upper insets. Parameter values are Vo = 7r/2, 
a = 2, mo = 1 and L = 50. The dashed areas are the stop 
bands of the corresponding Dirac-Kronig-Penney infinite lat- 
tice. 

the second grating region (at z > 0) comprising a se- 
quence of equally-spaced phase slips, at a distance a, 
equal to A<j) — 2Vo- The existence of surface states can 
be simply recognized by the appearance of narrow res- 
onance peaks embedded in a stop band region of the 
transmission spectrum of the grating (for more details 
we refer to [4"1"]V 

7 Photonic realizations of non-Hermitian 
relativistic wave equations 

Since the the seminal paper by Bender and Boettcher 
[51] . a great attention has been devoted toward the in- 
vestigation of non-Hermitian extensions of quantum me- 
chanics and quantum field theories. Indeed, many works 
have remarked that the Hermiticity of the underlying 
Hamiltonian can be relaxed, and that a a consistent 
quantum theory can be constructed for a broader class of 
Hamiltonians 85,86,87,88 , in particular those possess- 
ing parity-time (VT) symmetry. Non-Hermitian Hamil- 
tonians are also very often found in reduced descriptions 
of open Hermitian quantum systems, with important ap- 
plications to atomic, molecular and condensed-matter 
physics [89) . Several recent works have shown that op- 
tical structures in media with a complex refractive in- 
dex can provide an experimentally accessible test bed 
to simulate in a purely classical setting non-Hermitian 
features rooted in the non-relativistic Schrodinger equa- 
tion with a complex potential [9TI^l97M3"ll94ll95lig6] . 



Recently, non-Hermitian extensions of relativistic wave 
equations 97,98,99 and non-Hermitian quantum field 
theories |100j have attracted an increasing interest as 
well, however their physical realizations remain mostly 
unexplored. In Ref.[44 , optical simulations of non-Hermitian 
wave equations have been proposed, which are based on 
light propagation in distributed-fecdback (DFB) optical 
structures with controlled gain and/or loss regions. 

Let us indicate by n(z) = no — Anh(z) cos(2nz/ A + 
29(z)) the effective index grating of the DFB structure, 
where no is the modal refractive index in absence of the 
grating, An <C no and A are the peak index change and 
the nominal period of the grating, respectively, and h(z), 
29(z) are the normalized amplitude and phase profiles, 
respectively, of the grating. For a pure index grating, 
h{z) is real- valued, whereas for a pure gain grating h(z) 
is purely imaginary; in the most general case h{z) can 
be taken to be complex- valued [101] . The periodic mod- 
ulation of the refractive index leads to Bragg scatter- 
ing between two counterpropagating waves at frequen- 
cies close to the Bragg frequency u>b — 7rc/(Ano). The 
linear space-dependent absorption coefficient of counter- 
propagating waves in the structure is indicated by cxo(z) 
(cvo > in lossy regions, ao < in gain regions). In 
a semiconductor DFB structure, gain and loss regions 
could be tailored by a judicious control of current injec- 
tion across the active layer [101] . Indicating by E(z, r) = 
ipi(z, t) exp[— iu}BT+ikBZ+iO(z)]+ip2(z, r) exp[— iujbt— 
iksz — i0{z)] + c.c. the electric field propagating in the 
DFB structure, where ks = ^ / A, the envelopes ip\ and 
ip2 of counterpropagating waves satisfy coupled-mode 
equations [101] . After introduction of the scaled space 
and time variables x — z/Z and t = t/T, with Z = 
2noA/ '(-7T An) and T = Z/v g , where v g ~ c/no is the 
group velocity of light at frequency lub, the envelopes 
ipi and ip2 satisfy the following Dirac-type equation in 
the Weyl representation [33] 

id t i(j = —icr z d x ip + o x m(x)^ + V(x)ip = Hi/j (29) 

with complex-valued mass m and vector potential V 
given by 

d9 

m(x) = h(x), V(x) = — - i-y(x), (30) 
ax 

where 7(2:) = Zao(x) is the dimensionless absorption co- 
efficient and <j x ^ z are the Pauli matrices. In a DFB struc- 
ture, both m(x) and V(x) have a limited support over 
a spatial length L (the grating region), i.e. m = V = 
for \x\ > L/2. Since H is non-Hermitian, its spectrum 
is generally complex-valued, however it could happen 
that, in spite of non- Hermiticity, the energy spectrum 
remains real-valued, like in VT symmetric Hamiltoni- 
ans in the unbroken symmetry phase. In addition, be- 
cause of the non-self-adjointness of H, spectral singu- 
larities could arise in the spectrum of H, which corre- 
spond to cither a zero-width resonance in the transmis- 
sion/reflection spectrum of the DFB (the threshold for 
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self-oscillation) 94"ll4"4"] or to a perfect absorption of ra- 
diation under a suitable two-port coherent excitation of 
the DFB structure [411102] . Examples of non-Hermitian 
Dirac Hamiltonians showing such a two different kinds 
of spectral singularities were discussed in Ref.[H] as- 
suming a pure index grating (i.e. h(x) real-valued) for 
either 'PT-invariant and 'PT'-non-invariant DFB struc- 
tures. Another interesting example of a non-Hermitian 
relativistic wave equation is obtained by considering a 
purely gain-grating (h purely imaginary) in the absence 
of absorption losses and chirp (V = 0). In this case, for 
an infinitely- long grating (L oo, m = imo constant, 
with mo real- valued) Eq.(29) represents a superluminal 
extension of the relativistic Dirac equation [103] describ- 
ing a freely-moving hypothetical tachyonic particle, cor- 
responding to a negative mass square at rest (see e.g. 
[104 ). A remarkable property of such an equation, in 
addition to enable a kind of superluminal propagation, 
is the possibility to make a DFB structure fully transpar- 
ent, as recently discussed in Ref. |105| . Tachyonic exten- 
sions of Dirac equations can be realized as well in two- 
dimensional honeycomb photonic lattices made of cou- 
pled waveguides with alternating gain and loss regions 
(a complex photonic graphene), as recently proposed in 
Ref.gS]. 

8 Conclusion and outlook 

In this article, a brief overview on the possibility of- 
fered by light transport in periodic photonic structures 
to simulate in a purely classical setting the optical ana- 
logues of a wide variety of quantum phenomena rooted 
in relativistic wave equations has been presented. Spa- 
tial or temporal light transport in engineered photonic 
lattices and Bragg grating structures can simulate the 
Zitterbewegung of a relativistic electron, Klein tunnel- 
ing, vacuum decay and pair-production, the Dirac os- 
cillator, the relativistic Kronig-Penney model, and cer- 
tain non-Hermitian extensions of the Dirac equations, 
including superluminal (tachyonic) wave equations. Fur- 
ther quantum-optical analogies are expected to be in- 
vestigated and to be implemented in experiments using 
waveguide lattices and passive/active Bragg grating and 
DFB structures. 
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